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FOR THE STAR COLLABORATION 

Dynamic fluctuations in the local density of non-identified hadron tracks recon- 
structed in the STAR TPC are studied using the discrete wavelet transform power 
spectrum technique which involves mixed event reference sample comparison. The 
two-dimensional event-by-event analysis is performed in pseudo-rapidity n and az- 
imuthal angle <f> in bins of transverse momentum px ■ HIJING simulations indicate 
that jets and mini-jets result in characteristic signals, visible already at soft px, 
when the dynamic texture analysis is applied. In this analysis, the discrepancy be- 
tween the experiment and the HIJING expectations for Au+Au at \/Snn = 200 
GcV is most prominent in the central collisions where we observe the long range 
fluctuations to be enhanced at low px, and suppressed above px = 0.6 GeV 



1. Introduction 

The on-going RHIC program, motivated by an interest in the bulk prop- 
erties of strongly interacting matter under extreme conditions, has already 
yielded a number of tantalizing results. Deconfincmcnt and chiral symme- 
try restoration 1 are expected to take place in collisions of ultra-relativistic 
nuclei. Because these phase transitions are multiparticle phenomena, a 
promising, albeit challenging, approach is the study of dynamics of large 
groups of final state particles. The dynamics shows itself in the correlations 
and fluctuations (texture) on a variety of distance scales in momentum 
space. 

The multi-resolution dynamic texture approach (applied for the first 
time 2 at SPS) uses discrete wavelet transform 3 (DWT) to extract such 
information. At the present stage, the information is extracted in a com- 
prehensive way, without any built-in assumptions or filters. Mixed events 
are used as a reference for comparison in search for dynamic effects. Event 
generators are used to "train intuition" in recognizing manifestations of 
familiar physics (such as elliptic flow or jets) in the analysis output, as 
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well as to quantify sensitivity to the effects yet unidentified, such as critical 
fluctuations or clustering of new phase at hadronization. 

2. The STAR experiment 

The STAR Time Projection Chamber 4 (TPC) is mounted inside a solenoidal 
magnet. It tracks charged particles within a large acceptance (|?7| < 1.3, 
< (/> < 2ir) and is well suited for event-by-event physics and in-depth 
studies of event structure. The data being reported are obtained during 
the second (^/Snn — 200 GeV) year of RHIC operation. The minimum 
bias trigger discriminates on a neutral spectator signal in the Zero Degree 
Calorimeters 5 . By adding a requirement of high charged multiplicity within 
\r)\ < 1 from the scintillating Central Trigger Barrel, one obtains the central 
trigger. Vertex reconstruction is based on the TPC tracking. Only high 
quality tracks found to pass within 3 cm of the event vertex are accepted 
for the texture analysis. 

3. Dynamic texture analysis procedure 

Discrete wavelets are a set of functions, each having a proper width, or scale, 
and a proper location so that the function differs from only within that 
width and around that location. The set of possible scales and locations 
is discrete. The DWT transforms the collision event in pseudo-rapidity r\ 
and azimuthal angle (f> into a set of two-dimensional functions. The basis 
functions are defined in the (77, </3) space and are orthogonal with respect 
to scale and location. We accumulate texture information by averaging the 
power spectra of many events. 

The simplest DWT basis is the Haar wavelet, built upon the scaling 
function g(x) = 1 for < x < 1 and otherwise. The function 

f(x) = {+1 for < x < -; -1 for - < x < 1; otherwise} (1) 

is the wavelet function. 

The experimental acceptance in r],<j), and pt {{\ti\ < 1, < <f> < 2w)) is 
partitioned into bins. The rj-(f> partitions are of equal size, whereas in p?, 
the binning is exponential when more than one px bin is used. In each bin, 
the number of reconstructed tracks satisfying the quality cuts is counted. 

The scaling function of the Haar basis in two dimensions (2D) G{4>, rj) = 
g{4>)g{r]) is just a bin's acceptance (modulo units). The wavelet functions 
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Figure 1. Haar wavelet basis in two dimensions. The three modes of directional sen- 
sitivity are: a) diagonal b) azimuthal c) pseudo-rapidity. For the finest scale used, the 
white rectangle drawn "on top" of the function in panel a) would correspond to the 
smallest acceptance bin (pixel). Every subsequent coarser scale is obtained by expand- 
ing the functions of the previous scale by a factor of 2 in both dimensions. (Reproduced 
from 2 ). 

F x (where the mode of directional sensitivity A can be azimuthal (j>, pseudo- 
rapidity T], or diagonal <prj) are 

F*> = ffflfb), F* = f{<l>)g{Tj), F> = g(4>)f(v)- (2) 

We set up a two dimensional (2D) wavelet basis: 

F^jfarj) = 2™F\2 m cf>-i,2 m r,-j), (3) 

where m is the integer scale fineness index, i and j index the positions of 
bin centers in 4> and rj. Then, F^ i j with integer m, i, and j are known 3 to 
form a complete orthonormal basis in the space of all measurable functions 
defined on the continuum of real numbers L 2 (R). We construct G m .i.j(0, rj) 
analogously to Eq.3. 

Fig. 1 shows the wavelet basis functions F in two dimensions. At first 
glance it might seem surprising that, unlike the ID case, both / and g enter 
the wavelet basis in 2D. Fig. 1 clarifies this: in order to fully encode an 
arbitrary shape of a measurable 2D function, one considers it as an addition 
of a change along <f) (f((f>)g(rj), panel (b)), a change along rj (g{4>)f(rj), 
panel (c)), and a saddle-point pattern (/(</>)/ (77), panel (a)), added with 
appropriate weight (positive, negative or zero), for a variety of scales. The 
finest scale available is limited by the two track resolution, and, due to 
the needs of event mixing, by the number of available events. The coarser 
scales correspond to successively rc-binning the track distribution. The 
analysis is best visualized by considering the scaling function G m jj(<f),ri) 
as binning the track distribution p(<p, rj) in bins i,j of fineness m, while the 
set of wavelet functions F^ li A(j), rj) (or, to be exact, the wavelet expansion 
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)) gives the difference distribution between the data 



binned with given coarseness and that with binning one step finer. We use 
WAILI 6 software to obtain the wavelet expansions. 

In two dimensions, it is informative to present the three modes of a 
power spectrum with different directions of sensitivity P^ v (m), P^(m), 
P r, {m) separately. We define the power spectrum as 



where the denominator gives the meaning of spectral density to the ob- 
servable. So defined, the P x (m) of a random white noise field is indepen- 
dent of to. However, for physical events one finds P x (m) to be dependent 
on m due to the presence of static texture features such as acceptance 
asymmetries and imperfections (albeit minor in STAR), and non- uniformity 
of the dN/ dr] shape. In order to extract the dynamic signal, we use 
P x {m)true — P x {m)mix where the latter denotes power spectrum obtained 
from the mixed events. The mixed events are composed of the (to cf>) pix- 
els of true events, so that a pixel is an acceptance element of the finest scale 
used in the analysis, and in no mixed event is there more than one pixel 
from any given true event. The minimum granularity used in the analysis 
is 16 x 16 pixels. a 

Systematic errors can be induced on P x (m)true — P X ('m)mix by the pro- 
cess of event mixing. For example, in events with different vertex position 
along the beam axis, same values of 77 may correspond to different parts of 
the TPC with different tracking efficiency. That will fake a dynamic tex- 
ture effect in to In order to minimize such errors, events are classified into 
event classes with similar multiplicity and vertex position. Event mix- 
ing is done and P x (m) trU e — P X (m) m i X is constructed within such classes. 
Only events with z vertex lying on the beam axis within 25 cm from the 
center of the chamber are accepted for analysis. To form event classes, 
this interval is further subdivided into five bins. We also avoid mixing of 
events with largely different multiplicity. Therefore, another dimension of 
the event class definition is that of the multiplicity of high quality tracks 
in the TPC. For central trigger events, the multiplicity range of an event 
class is typically 50. 




(4) 



a For a quick reference, here are the scales in r). Scale 1: Arj = 1; scale 2: At] = 1/2; 
scale 3: Ar) = 1/4 and so on. 
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4. "Coherent" interference of patterns and normalization 
of power spectra 

Imagine a reconstructed event as a distribution of points in the space of 
variables (j],<f>,px)- We slice this space into px bins and analyze two- 
dimensional (77, 4>) patterns. The patterns from different px slices of the 
same event will amplify the texture signal when those px bins are merged. 
Depending on how the amplification works, one will find different scaling 
law to relate the P A (m) true — P x {m) m i X signal amplitude with the under- 
lying number of particles. 

The DWT power spectrum at each scale is (using Haar wavelet) a sum of 
squared pixel-to-pixel content differences for the given pixel fineness (scale) . 
One can think of the pixel-to-pixel content difference the same way as one 
thinks of a random fluctuation in the pixel content. Imagine that the pat- 
tern being analyzed is a small sub-sample of the event, and its number of 
particles N can be increased at will, up to the point of making it an en- 
tire event - as is the case when the sub-sample is a px bin of the event. 
The pixel content will scale with N, and if the dynamic pattern preserves 
its shape from one px bin to another, the pixel-to-pixel difference on the 
characteristic scale of the pattern will also scale as N. Consequently, the 
dynamic component of the power spectrum for this scale will grow as N 2 . 
We will call this behavior "coherent" in analogy with optics, where one 
needs coherence in order to see interference patterns. 

Normalization is needed in order to, first, express different measure- 
ments in the same units; second, eliminate trends in px dependence which 
are induced by the design of the measure and unrelated to the physics. For 
the "coherent" case, the normalized dynamic texture observable is 

(P X (m)true ~ P X {m) mix )/P x {m) mix /N. (5) 

One could also imagine "incoherent" px slices. In the "incoherent" case, the 
pixel content will grow proportionally to N, but the pixel-to-pixel difference 
will grow as the RMS fluctuation of the pixel content, i.e. as the Poissonian 
s/~N. The dynamic component of the power spectrum will grow as N (i.e. 
oc P{m)) and 

{P X {m)t r ue ~ P X (m) mix )/P X (m) mix (6) 

should be used in this case. In the DWT-bascd texture analysis, amplifi- 
cation of the signal is based not on adding the patterns themselves, but on 
adding the power spectra of local density fluctuations, that is (continuing 
the optics analogy) adding the intensities rather than field amplitudes. For 
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Figure 2. a) - rj(j>, b) - <j> and c) — ij 
directional components of the dynamic 
texture in HIJING (events with impact 
parameter between and 3 fm), arising 
primarily due to jets. Data sets with dif- 
ferent pt bin widths, indicated by the 
open and solid symbols, are statistically 
consistent at both scales when the "co- 
herent" normalization is included. • o— 
scale 1; • o — scale 2. Enhanced fineness 
scale 2 of the <p texture plot (b) reflects 
back- to-back correlations. 
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this reason, in the DWT analysis one does not require "coherence" to am- 
plify the signals from many px slices, just as in optics one does not need 
coherence to see the light intensity increase with an increase in the number 
of photons. 

5. Textures of jets and critical fluctuations in event 
generators 

Dynamic texture is to be expected from HIJING 7 given its particle produc- 
tion mechanism at RHIC energy (jets, mini-jets and string fragmentation). 
HIJING combines a perturbative QCD description at high pr with a model 
of low pt processes. Figure 2 demonstrates observability of the HIJING 
dynamic effects in our analysis. b We see that, first, the difference between 
the true and mixed events is noticeable and can be studied as a function of 
Pt with the present HIJING statistics of around 1.6 x 10 5 events. Second, 



b In all MC generators, no GEANT and no response simulation is done. Instead, only 
stable charged particles (e,fi,w,K,p) and their antiparticles from the generator output 
are considered, provided that they fit into the STAR TPC fiducial r\ acceptance |tj| < 1. 
Momentum resolution and pt acceptance are not simulated. 
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Figure 3. (P£ ue - 

PliJ/PliJ* from the Crit- 
ical MC generator. Events 
with 20 to 30 charged tracks in 
the STAR acceptance are ana- 
0.1 0.2 0.3 0.4 0.5 lyzcd. O ~~ scale 1, - scale 2, 

p T- GeV/0 A - scale 3, □ - scale 4. 

the open and closed symbols, which correspond to different px bin sizes, ap- 
pear to fall on the same curve after the 1/N normalization, where N is a px 
bin multiplicity, as would be the case for "coherent" (see Section4) px bins. 
Third, the rise of the signal with pt is due to the fact that high px is dom- 
inated by jet production. As far as the px "coherence" is concerned, one 
would expect that a high px parton, creating hadrons via fragmentation, 
produces similar {q,<j)) patterns at different px as the energy sharing among 
the secondaries proceeds, and thus the coherent interference of px patterns 
is natural for this mechanism of particle production. These signals in HI- 
JING arc gone when jet production is turned off in the generator. Ability to 
study jet textures at soft, as well as high, px means that the study promises 
to be very informative because majority of the reconstructed tracks will be 
utilized. 

CMC is Critical Monte Carlo generator created by N.Antoniou and 
coworkers 8 . In the framework of an effective action approach, these authors 
simulate a system undergoing a second order QCD phase transition. The 
r\ signal at low px (Fig. 3) is much stronger than seen in HIJING and is 
dominated by the coarse scale. 

6. STAR measurements of dynamic textures 

Elliptic flow is a prominent large scale dynamic texture effect already well 
measured at RHIC 9 . The DWT approach localizes elliptic flow on scales 
2 and, to some degree, 3 of the azimuthal observables. In this report, we 
ignore flow and concentrate on the r\ observables. 

Fig. 4 presents the STAR measurements of long range (scale 1) fluc- 
tuations in peripheral (0.014 < mult/no < 0.1) collisions and compares 
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Figure 4. (Pt rue — ^mix) / ^mix / ^ f° r sca -l e b peripheral events. Open stars - STAR 
data for ^fs = 200 GeV, 0.014 < mult/n < 0.1. • - HIJING at the same energy, 
mult/no < 0.1. 



them with HIJING simulations. Qualitatively, both sets of points behave 
similarly: a region of nearly flat or falling behavior around mean px is 
replaced by a rising trend for px > 0.8 GeV/c. This trend has already 
been discussed in Section 5 and is due to jets. The HIJING signal is below 
the STAR data at low pr, but reaches higher values at higher px\ its rise 
with pt is stronger. From this figure we conclude that the fluctuations in 
local hadron density due to jet production are observable at RHIC in the 
soft pt range (p? < 2 GeV), and that their qualitative features are rea- 
sonably well described by a super-position of independent nuclcon-nucleon 
collisions based on the physics learned from pp(p) and e + e~ experiments at 
comparable energies. Quantitatively speaking, we keep in mind that due to 
nuclear shadowing effect 11 , peripheral Au+Au events are not supposed to 
be identical to elementary collisions. A comparison of pp, dAu and AuAu 
data from RHIC will shed more light on this effect. In the absence of exper- 
imental data on nuclear shadowing of gluons, HIJING assumes 7 equivalence 
of the effect for quarks and gluons. 

Next look at a central sample (Fig. 5) - there is a remarkable differ- 
ence: we now see a change in the pr trend above pr — 0.6 GeV. Instead of 
rising with pr (as in the peripheral events), the STAR data points become 
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Figure 5. (P? rue - Pmix)l P mixl N for scale central events (0.65 < mult/no < 1.). 
Open stars - STAR data for = 200 GeV. • - regular HIJING; o - HIJING with jet 
quenching, both at ^/s = 130GeV. 



consistent with 0. The pr trends in the data and HIJING look opposite: 
the model still predicts a monotonic rise with px- Can there be a single 
explanation to both disappearance of texture at moderate pt and its en- 
hancement at low pt? The hypothetical deconfined medium is expected to 
suppress jet production via dissipative processes (jet quenching) 10 . The 
medium-induced energy loss per unit of length is proportional to the size of 
the medium and thus, the effect grows non-lincarly with system size. Sup- 
pression of hadron yields at high px in central AuAu events with respect to 
scaled pp and peripheral collisions has been reported 12 and interpreted as 
an evidence of medium effects (possibly, nuclear shadowing 11 ). Jet quench- 
ing is modeled in HIJING, and is seen (compare two sets of HIJING points 
in Fig. 5) to affect the texture observable somewhat. If the dissipation takes 
place, one may expect that as jets and mini-jets thermalize, the textures 
associated with them migrate towards mean pr- A transport model would 
be needed in order to simulate such a process. However, the low pT fluctu- 
ations may have an independent origin, unrelated directly to the partonic 
energy loss in medium. 
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7. Conclusions 

A non-trivial picture of texture effects emerges when the DWT power spec- 
trum technique is applied to AuAu data from RHIC. Long range (Ar/ wl) 
pseudo-rapidity fluctuations at soft px are observed in peripheral events 
and identified with jets and mini-jets. In central events, these fluctuations 
are not seen, which indicates a change in the properties of the medium. 
Large scale of the effect points to its early origin. An excess of fluctuations 
at low pt compared to HIJ1NG is seen in peripheral and central events. 
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